Dissipative processes in superfluid quark matter 
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Abstract. We present some results about dissipative processes in fermionic superfluids that are 
relevant for compact stars. At sufficiently low temperatures the transport properties of a superfluid 
• are dominated by phonons. We report the values of the bulk viscosity, shear viscosity and thermal 

I conductivity of phonons in quark matter at extremely high density and low temperature. Then, we 

present a new dissipative mechanism that can operate in compact stars and that is named "rocket 
term". The effect of this dissipative mechanism on superfluid r-mode oscillations is sketched. 
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INTRODUCTION 

Superfluidity is a property of quantum fluids related with the existence of low energy 



excitations that satisfy the Landau's criterion for superfluidity [1, 2, 3] 
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where e(p) is the dispersion law of the excitations. The main property of superfluids is 
that they have an extremely small viscosity and can flow with very small friction. 

In general, superfluidity is due to the appearance of a condensate which spontaneously 
breaks a global symmetry of the system. As a consequence, in the low energy spectrum 
one has a Nambu-Goldstone boson, <p, whose dispersion law satisfies Eq. (1). 

Superfluidity was first discovered in 4 He, which becomes superfluid when cooled at 
temperatures below 2. 17 K. The superfluid property of 4 He atoms is due to the fact that 
they are bosonic with very small interactions. Therefore at low enough temperature they 
Bose-Einstein condense in the lowest quantum state. 

Fermionic systems can become superfluid as well. According with the Cooper theo- 
rem, fermionic superfluidity takes place in quantum degenerate systems when the inter- 
action between neutral fermions is attractive and the temperature is sufficiently low. In 
this case one has the formation of di-fermion condensates, or Cooper pairs, that behave 
as bosonic particle and can condense. Also in this case superfluidity is accompanied 
with the spontaneous breaking of a global U(l) symmetry and with the appearance of a 
Goldstone boson. 

Fermionic superfluids might be present in the interior of neutron stars where the tem- 
perature is low and the typical energy scale of particles is extremely high. In particular, 



in the inner crust of neutron stars the attractive interaction between neutrons can lead to 
the formation of a BCS condensate. Moreover, if deconfined quark matter is present in 
the core of neutron stars it is very likely to be in a color superconducting phase [4]. 

A means to detect or to discard the presence of a superfluid phase in compact stars 
consists in studying the evolution of the r-mode oscillations [5]. R-modes are non-radial 
oscillations of the star with the Coriolis force acting as the restoring force. They provide 
a severe limitation on the rotation frequency of the star through coupling to gravitational 
radiation. When dissipative phenomena damp these r-modes the star can rotate without 
losing angular momentum to gravitational radiation. If dissipative phenomena are not 
strong enough, these oscillations will grow exponentially and the star will keep slowing 
down until some dissipation mechanism is able to damp the r-mode oscillations. In this 
way one can put some constrains on the stellar structure ruling out phases that do not 
have large enough viscosity. For such studies it is necessary to consider in detail all the 
dissipative processes. 



The transport properties of a superfluid depend on the shear viscosity coefficient, rj, on 
three independent bulk viscosity coefficients, £1, £2, £3, and on the thermal conductivity 
K. Microscopically, these coefficients are related with interactions that take place in the 
superfluid. Macroscopically, they can be viewed as the phenomenological coefficients 
which relate the rate of change of currents with the corresponding affinities [6]. The 
requirement that the dissipative processes lead to positive entropy production imposes 
some constraint on these coefficients: K, 7], £2, £3 must be positive and £ 2 < £2 £3- 

The presence of three different bulk viscosity coefficients might sound puzzling. It is 
related with the fact that the hydrodynamic equations governing a superfluid are essen- 
tially different from standard fluid equations. In a superfluid there are two independent 
motions, one normal and the other superfluid. While the coefficient £2 plays the role 
of the standard bulk viscosity coefficient, the coefficients £1 and £3 provide a coupling 
between the normal and superfluid components. The friction forces due to bulk viscosi- 
ties can be understood as drops, with respect to their equilibrium values, in the main 
driving forces acting on the normal and superfluid components. These forces are given 
by the gradients of the pressure P and of the chemical potential /1. One can write in the 
comoving frame 



where P eq and jU eq are the equilibrium pressure and chemical potential, V is a quantity 
proportional to the quantum condensate, co^ = — (d^(p + jiu^ 1 ) and uP is the velocity of 
the fluid. 

In a conformally invariant system it has been shown in Ref. [7] that £1 = £2 = 0. 
However, £3, k and r\ cannot be determined by the same symmetry argument. 

In the low temperature regime, T <C T c , where T c is the critical temperature for 
superfluidity, the transport properties of superfluids are determined by phonons. The 



DISSIPATIVE PROCESSES IN SUPERFLUIDS 



P = J P eq -Cidiv(y 2 w)-C 2 divu, 
H = i Ueq-C3div(V 2 w) -Cidivu, 



(2) 
(3) 



contribution of other degrees of freedom is thermally suppressed. In this case one 
can show that Q = ^£3, meaning that there are only two independent bulk viscosity 
coefficients. 

For T <C T c the transport coefficients strongly depend on the phonon dispersion law. 
The shear viscosity is the only transport coefficient that does not vanish for phonons with 
a linear dispersion law. But, for the bulk viscosities and for the thermal conductivity one 
has to include the term cubic in momentum 

e(p) = c s p + Bp 3 + 0(p 5 ). (4) 

Moreover, in the computation of the bulk viscosity one has to consider the processes that 
change the number of phonons. The parameter B determines whether some processes 
are or are not kinematically allowed. For B > the leading contribution comes from the 
Beliaev process <p — > (p(p. In the opposite case the Beliaev process is not kinematically 
allowed and one has to consider the processes <p<p — > (pqxp. 

A different dissipative process is associated with the change in the number densities 
of the various species present in the superfluid, which leads to the appearance of the so- 
called rocket term in the Euler equations of the system. This force is due to the fact that 
when two or more fluids move with different velocities a change of one component into 
the other results in a momentum transfer between the fluids. This change in momentum 
is not reversible, because it is always the faster moving fluid that will lose momentum. 
The resulting dissipative force is proportional to the mass rate change, and to the relative 
velocity of the fluids. The name "rocket term" reminds that the same phenomenon takes 
place in the dynamical evolution of a rocket whose mass is changing in time as it 
consumes its fuel. In compact stars number changing processes can take place in the 
outer core and in the inner crust of the star and are related to the Urea processes and/or 
interactions between the superfluid and the crust. The dissipative force due to the rocket 
term has been considered in the context of r-mode oscillations of standard npe neutron 
stars, for the first time in Ref. [8]. 



COLD AND DENSE QUARK MATTER 

Quantum Chromodynamics predicts that at asymptotically high densities quark matter 
is in the color- flavor locked (CFL) phase [9]. In this phase up, down and strange quarks 
of all three colors pair forming a di-fermion condensate that breaks the U ( 1 )b symmetry 
spontaneously. 

At asymptotic densities, the dispersion law of phonons in the CFL phase has been 
derived in [10], and one has that c s = and B = — , where A is the gap parameter. 
The effect of the quark masses can be taken into account perturbatively, as far as m s <C /i, 
moreover the coupling constant is small, <C 1, and one can assume that CFL is 
approximately scale invariant. 

A first study of the shear viscosity and of the contribution to the bulk viscosity 
coefficients due to phonons has been presented in Refs. [11, 12, 13]. Beside phonons, 
kaons may give a sizable contribution to the transport coefficients. The contribution of 



kaons to £2 has been studied in Ref. [14]. There is still no computation of the contribution 
of kaons to the remaining bulk viscosity coefficients. 

Considering only the contribution of phonons, in Ref. [13] it has been found that 
neglecting the running of the QCD gauge coupling 
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Notice that for vanishing quark masses the system is scale invariant and £1 and £2 vanish. 

The contribution to thermal conductivity due to phonons and kaons has recently been 
studied in Ref. [15]. The thermal conductivity from phonons turns out to be dominant 
and given by 

K~6xl(T 2 ^, (6) 
whereas the shear viscosity from phonons as determined in [1 1] is given by 
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T^OxKT 2 ^. (7) 

A different dissipative mechanism is due to mutual friction, that has been studied in 
the CFL phase in Refs.[16] and [17]. 



ROCKET TERM 

For a system consisting of neutrons, protons and electrons, the mass conservation law is 
given by, see e.g. [18], 

d t p x +v l (p x v i x )=r x , (8) 

where T x is the particle mass creation rate per unit volume and the index x = n,p,e refer 
to the particle species, that is, neutrons, protons and electrons. In these equations we 
have considered that some process can convert neutrons in protons and electrons and 
vice versa. Therefore, we are assuming that the various components are not separately 
conserved. One possible mechanism leading to a change in the particle number densities 
is given by the Urea process 

«->■ p + e~ + v e , p + e~ —>n + v e . (9) 

A different process is the crust-core transfusion, a process that takes place when the 
ionic constituent of the crust are squeezed by the underlying superfluid and part of their 
content is released and augments the superfluid component. The opposite mechanism, 
related to a reduction of the pressure leads to the neutron capture by the ions of the crust. 
We restrict our analysis to the beta decay processes. 

The three particle creation rates are not independent quantities, because charge and 
baryon number conservations imply that T e = T p — T n , meaning that only one creation 
rate is independent. 

It is possible to simplify the treatment of the system considering that electrons and 
protons are locked together by the electromagnetic interaction. Thus they move with the 
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FIGURE 1. Instability window of the superfluid r-modes of a star with uniform density p = 2.5po, with 
radius R = 10 km and mass M ~ 1.47M©. The dashed red line represents the instability window in the 
absence of the rocket term. The full blue line represents the instability window with the inclusion of the 
rocket term. The horizontal full lines correspond to the effect of the mutual friction for £ = 0.0002, lower 
line, and e = 0.002 upper line. In our simplified model of star, the mutual friction is independent of the 
temperature. 



same velocity. Moreover, charge neutrality implies that n e = n p . Therefore, electrons 
and protons can be treated as a single charge-neutral fluid and henceforth we shall refer 
to this fluid as the "charged" component. The two fluids have mass densities 

p n = m n n n and p c = m n n c , (1 0) 

where m n = m p + m e and n e = n p = n c . The Euler equations obeyed by the two fluids 
are given by 

(4 + viV i )(v? + £ n w / )+V i (/i )I + $)+^ = 0, (11) 
(^ + v^'V y )(v?-e c w l -)+V l -(/i c + ^)-e c ^'V / v5 = (1 - e n - e c ) ^ , (12) 

Pc 

where <J» is the gravitational potential, i, j label the space components; we have defined a 
chemical potential by mass \x x = ll x /m n , and w = v c — v„ represents the relative velocity 
between the two fluids. The quantities e n and e c are the entrainment parameters, that are 
related to the fact that momenta and velocities of quasiparticles may not be aligned [19]. 

The term on the right hand side of Eq. (12) is the rocket term. In the analysis of 
the possible dissipative mechanisms of star oscillations this term is usually neglected. 
Indeed, it is in general assumed that the neutron, proton and electron numbers are 
separately conserved quantities, that is T p = Y e = Y n = 0. In Ref. [8] we have considered 
the effect of this term and obtained the corresponding damping timescale, assuming that 
the change in the number densities is due to direct Urea processes. 

Stability window for superfluid r-modes 

In Fig. 1 we report the result for the superfluid r-mode "instability window" for a star 
with uniform density, p = 2.5po, and R = 10 km. It is obtained comparing the various 



dissipative timescales with the growth timescale of gravitational wave radiation. Stars 
above the various lines are unstable. The dashed red line represents the instability win- 
dow in the absence of the rocket term. The region above the dashed red line is unstable 
when only shear and bulk viscosity damping are considered. At low temperature, shear 
viscosity is the dominant dissipative mechanism. With increasing temperature shear vis- 
cosity is less efficient and it can damp r-mode oscillations for smaller and smaller values 
of the frequency. The behavior of bulk viscosity is the opposite and starts to damp r- 
mode oscillations for temperatures of the order of 10 10 K. The full blue line represents 
the effect of the rocket term. It has a behavior qualitatively similar to the bulk viscosity, 
but it becomes effective at smaller temperatures. Therefore, the instability window with 
the inclusion of the rocket term is much reduced. The full horizontal lines correspond 
to the effect of mutual friction, see e.g. [20], for two different values of the entrainment 
parameter. 
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